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We propose a simple dynamic model of polymers under shear with an anisotropic mobility tensor. 
We calculate the shear viscosity, the rheo-dielectric response function, and the parallel relaxation 
modulus under shear flow deduced from our model. We utilize recently developed linear response 
theories for nonequilibrium systems to calculate linear response functions. Our results are qual- 
itatively consistent with experimental results. We show that our anisotropic mobility model can 
reproduce essential dynamical nature of polymers under shear qualitatively. We compare our model 
with other models or theories such as the convective constraint release model or nonequilibrium 
linear response theories. 
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I. INTRODUCTION 

The linear response theory gives a formula which relates a response function in equilibrium to an equilibrium 
time correlation function [lH3|. I n equilibrium, the response of a physical quantity A at time t to a weak external 
perturbation at time t'(< t), which is conjugate to B, is given as the following form. 

K AB (t-t') = 1 ± f ^.{A(t)B(t')) eq (1) 

Here ks is the Boltmzann constant, T is the absolute temperature, and (• • • ) oq means the equilibrium statistical 
average. Eq ((T|) holds for a wide range of systems [H-IH as long as the system is in equilibrium. (In the followings, 
we call the response formula (p} as the Green-Kubo type formula.) From the view point of experiments, eq ([TJ can 
be utilized to obtain the correlation function from the response function. This enables us to extract information of 
microscopic and/or mesoscopic dynamics of molecules from macroscopic responses. 

For polymeric systems, the viscoelastic response functions (especially the shear relaxation modulus or the storage 
and loss moduli) or the dielectric response function are useful. For example, the viscoelastic response functions can 
be related to the autocorrelation function of the microscopic stress tensor, and the stress tensor can be related to 
bond vectors [4]. The dielectric response functions of polymers with type A dipoles (polymer chains which have 
electric dipoles along the chain backbones) can be related to the autocorrelation functions of end-to-end vectors. The 
combination of the viscoelastic and dielectric measurements provides various information about dynamics of polymer 
chains 

Even out of equilibrium, it is possible to measure linear response functions. Then we expect that the measured 
response functions can be related to the correlation functions, or at least they reflect the information about the 
dynamics of polymer chains. Actually, several linear response measurements under steady shear, such as mechanical 
responses flOl4l7j or dielectric responses [l8l - l22T ] have been reported. However, unlike the equilibrium cases, it is not 
clear how we can analyze and interpret the nonequilibrium linear response functions. In several works, the Green- 
Kubo type relation ([1} is utilized to analyze obtained experimental data, without justifications. But from the view 
point of nonequilibrium statistical physics, generally eq ^ does not hold except for some special or limited cases. 

Fortunately, the Green-Kubo type relation approximately holds for the dielectric response of polymers in the shear 
gradient direction under shear (the rheo-dielectric response) (23|, and thus we can obtain the correlation function of 
the end-to-end vectors under shear. In this work we therefore concentrate on polymers with type A dipoles. The 
rheo-dielectric responses for polymers with type A dipoles has been systematically studied and analyzed for various 



systems and shear rates [18j, l2Cj-|22j , The rheo-dielectric functions of linear polymers are reported to be insensitive to 
shear, even if the shear rate is large and the shear thinning is observed. This implies that the dynamics of polymer 
chains in the shear gradient direction is not affected largely by shear flow even under fast shear. So far, why and 
how such insensitivity occurs is not fully understood. Although attempts have been done by coarse-grained molecular 
simulations [23, 24], simulations have failed to reproduce experimentally obtained rheo-dielectric behavior. 

In the field of the constitutive equation models, many different approaches to nonequilibrium systems have been 
utilized, to reproduce nonlinear viscoelasticity well. Among them, the anisotropic mobility type models [25| - [30j are 
particularly notable. In the anisotropic friction model, the friction coefficient (or almost equivalently, the relaxation 
time) of the model is expressed as a tensor quantity instead of a scalar quantity. This allows us to reproduce variety 
of models which can reproduce complex viscoelastic behavior with relatively simple constitutive equations. 

Motivated by the anisotropic mobility type constitutive equation models, in this work we aim to propose a Langevin 
equation model with an anisotropic mobility tensor for dynamics of polymers under shear. Although the anisotropic 
mobility tensors are not utilized widely in the field of the nonequilibrium statistical mechanics, they can potentially 
model the dynamics under fast shear. A Langevin equation model with an anisotropic mobility tensor model was 
already proposed for a rather simple system [31[ based on the projection operator method and molecular dynamics 
simulation data. We construct an anisotropic mobility model which is consistent with some molecular dynamics simu- 
lation results, in a similar way. We limit ourselves to an analytically solvable toy model to examine the characteristic 
properties of the anisotropic mobility model explicitly and simply. We explicitly calculate several linear response 
functions such as the rheo-dielectric response function and the moduli in the shear flow direction (parallel moduli) 
for our model. To calculate the linear response functions, we utilize recently developed linear response theories for 
nonequilibrium Langevin systems (32T - [34j | . Finally we discuss the properties of our model and its linear response 
functions under shear from several aspects. We compare our model with several pieces of previous work, and show 
the differences and similarities between our model and the other models. 
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II. MODEL 



In this section, we propose a simple and solvable model for dynamics of polymers under shear. We consider weakly 
entangled polymer melts or solutions, for which rheo-dielectric experiments have been carried out fTsL [20I - I22I ] . There 
are several different approaches to model the dynamics of polymers 0, [H, H(| . Thus at first, we should choose an 
appropriate model from candidates. We are interested in the qualitative and essential feature of the polymer dynamics 
under shear. We require the model to be simple and analytically solvable, so that we avoid unnecessary complexities 
and confusions, and make the model properties clear. From these requirements, we limit ourselves to the dynamics 
of the end-to-end vector of a single polymer chain R, in the weakly entangled system. This reduces the degrees of 
freedom drastically. (This approximation gives so-called the dumbbell model [3y| .) We also limit ourselves to the 
dynamics in the long time limit, where is no memory effect for the end-to-end vector dynamics. Then we can describe 
the dynamic equation for R in a closed, Markovian form. In equilibrium, this equation can be described as 

dR(t) I dT(R(t)) 

where £0 is the friction coefficient which the end-to-end vector feels, J-{R) is the free energy, and £o(t) is the Gaussian 
noise. For a Gaussian chain, the free energy can be expressed in the following simple linear elasticity form. 

F(R) = ^R* (3) 

R 2 is the equilibrium mean square average end-to-end distance. The fluctuation-dissipation relation of the second 
kind [J requires for the thermal noise £o(t) to satisfy the following relations. 

<€<>(*)> = (4) 

(umo(t')) = ^s( t -t')i ( 5 ) 

Here (• • ■ } denotes the statistical average and 1 is the unit tensor. 

We consider a system under simple shear. The velocity gradient tensor n is given as follows. 

k -I** ( a = x >P = y) 

I (otherwise) 

where 7 is the shear rate. We assume that 7 is smaller than r^ 1 (tr is the Rouse time of the polymer chain) and the 
polymer chain is not so highly stretched. A schematic image of the system is shown in FIG. [T] Since the shear flow 
cannot be expressed in relation to a conservative potential force, the system is nonequilibrium under simple shear. 
As a simple extension of the equilibrium Langevin equation @, one may consider the following Langevin equation. 

— - - Co ~^wr + ■ m + m (7) 

In the Langevin equation ([7|l. the mobility, which is defined as the inverse of the friction coefficient, is scalar and thus 
isotropic. This is because we simply used the mobility (the friction tensor) in equilibrium. The system is isotropic 
in equilibrium, and the mobility should also be isotropic from the symmetry. However, if the system is not in or 
near equilibrium, it is generally not isotropic. This means that, in prin ciple, the mobility could be an anisotropic, 
tensor quantity [3l[ . Actually, in some constitutive equation models [27], Ho[ , the mobility tensors are designed to be 
anisotropic under flow. Thus we consider that eq (JTJ) is an oversimplified dynamic equation model, which may lead 
physically incorrect results in particular under fast flow. 

When the mobility tensor becomes anisotropic, we expect that some transport coefficient tensors also become 
anisotropic [37- 39] . Anisotropic diffusion coefficient tensors are actually observed in nonequilibrium molecular dy- 
namics (NEMD) simulations based on the SLLOD model [40l l4l[ . The diffusion tensor of Lennard- Jones (LJ) particles 
under shear is first reported by Sarman, Evans, and Baranyai [401 ] . They studied the diffusion tensor of an LJ particle 
under shear, at various rates. Their simulation result clearly shows that the diffusion tensor becomes anisotropic 
and it depends on the shear rate. Quite recently. Hunt and Todd 41] performed NEMD simulations for relatively 
short polymer melts (the Kremer-Grest chains [421 ] ) and measured the diffusion tensors under shear. Their result is 
qualitatively similar to the case of the LJ system. Namely, the diffusion tensor of the center of mass of a polymer 
chain becomes anisotropic and depends on the shear rate. Besides, the dependence of the diffusion tensor on the shear 
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rate becomes strong as the polymerization index (number of beads in a chain) increases. Although their simulations 
are limited for rather short chains (unentangled chains), we expect that this trend will be qualitatively the same or 
even enhanced for well entangled chains. 

Although there are several possible interpretations for these NEMD results [3l],[43j], in this work we interpret that the 
anisotropic diffusion is caused by anisotropic mobilities. Namely, if we model the dynamics of a polymer chain under 
shear by a Langevin equation, we should employ an anisotropic mobility tensor (or an anisotropic friction tensor). 
McPhie et al |31| proposed a Langevin equation with an anisotropic friction tensor to describe the coarse-grained 
motion of an LJ particle under shear. Although they proposed an underdamped Langevin equation, an overdamped 
Langevin equation (like eq ([2)1) seems to be more suitable for the end-to-end vector of a polymer. Then we express 
the Langevin equation under shear as follows. 

«._ A(7) .2£ffl + K . fl + « 7;t) 

where A(7) is the mobility tensor which depends on the shear rate 7 and £(7, t) is the Gaussian noise. We assume 
that the fluctuation-dissipation type relation between £(%£) and A(7) is satisfied. 

<€(75*)> = (9) 
(€(7; *)€(75 0> = 2fc B TA( 7 )<5(t - (10) 

Eq ([TO]) requires the mobility tensor A( 7 ) to be symmetric under the shear field given by eq ©■ Quite recently, Ilg 
and Kroger [44| IH| proposed a constitutive equation model with anisotropic mobility (friction coefficient) tensor for 
relatively short polymer chains, based on the NEMD results. Although their model is not equivalent to ours, it is 
qualitatively similar. 

From NEMD simulation results and properties of a Langevin type equation, we consider the mobility tensor A (7) 
under steady shear should have the following properties. 

• The mobility tensor should be positive definite (all of its eigenvalues should be positive). 



The eigenvalues of A( 7 ) are unchanged under the transform 7 — > —7. That is, each eigenvalue is an even 
function of 7. 

The xx-component (A xx ) decreases as 7 increases, while other diagonal components (A^ and A zz ) are not so 
sensitive to 7. 

The xy-component (A xy ) is non-zero but its value is smaller than diagonal elements. Thus it may be simply 
neglected (A xy « 0). 

From the symmetry, A xz — A yz = 0. 

The mobility tensor should reduce to the isotropic tensor at equilibrium (7 = 0). 



In this work, we employ the following simple form which satisfies the above properties. (We discuss about the other 
possible forms later.) 



A(7) = r 



A( 7 ) 
1 
1 



(11) 



Here, Co is the friction coefficient at equilibrium and A( 7 ) is a function of 7. A( 7 ) is an even function of 7, and it 
monotonically increases as j 2 increases. Further, since it should recover the equilibrium form in the absence of shear 
flow, A( 7 ) reduces to the equilibrium form, A(0) = 1, at the limit of 7 — > 0. For example, we can employ the following 
simple form for A (7). 

A(7)=[l + (rc7)T /2 (12) 

a is an exponent which satisfies < a < 1, and r c is a characteristic crossover time. Roughly speaking, the mobility 
tensor is isotropic for |r c7 | <C 1 and anisotropic for |t c 7| 3> 1. As we show in the next section, eq (TT2")) gives the 
power law type behavior of the shear viscosity. We call eq (|12[) as the power law type model. We can also employ the 
following form. 



(13) 
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Coo corresponds to the effective friction coefficient in the x-direction for |t c 'V| 3> 1. Eq (fl3l) may be preferred if the 
system exhibits the second Newtonian region for the shear viscosity. Eqs 1)120 and (TTBI are just possible candidates, and 
there are many other possible forms for A (7 ). It is worth noting that our anisotropic mobility tensor is qualitatively 
similar to the model by Ilg and Kroger |45| (in their model, the dynamics of the x-direction is also accelerated 
under shear). Before we proceed, we should notice that the anisotropic mobility tensor (fTTj) is designed for the simple 
shear flow, and it is not expected to be applicable for other flows such as elongational flows. In the following analysis, 
we consider the response around the steady state under the simple shear (expressed by eq (J6))) and thus eq (jlip is 
sufficient for our purpose. 

Since the Langevin equation (|8j) is linear in R, we can analytically integrate it and obtain explicit expression for 
several physical quantities. Thus we can analyze linear responses such as the rheo-dielectric response of our model 
explicitly. By substituting eq (fTTj) into eq (jSJ, we have the following equations for R x ,R y , and R z . 



TQ 



1 



dRxjt) = Mj) 
dt 

dRyjt) 

dt 
dR z {t) 
dt 



R x (t)+jR y {t) + Ui,t) 



= — Ry(t) + Ui,t) 

TO 

= --R z (t) + Ui,t) 

TO 



(14) 
(15) 
(16) 



where we defined the characteristic time tq. 



tq 



CoR 2 



3k B T 

The probability distribution defined via the following equation is useful for some calculations. 

P(r,t) = (6{r-Wt))) 



(17) 



(18) 



The probability density (fTBf follows the Fokker-Planck equation, which describes the time evolution due to the 
Langevin equation (|8]l. 



dt Or { 11 



9 ^P(r,t) + k B T dP ^ 



dr ' 1 dr 
The steady state solution of the Fokker-Planck equation (fT5|) . P ss (r) satisfies 



d 

-Q^-[K-rP(r,t)} 



(19) 



-j-P^ + ksT-^ 



d 

— ■ [k ■ rP ss {r)} 



(20) 



For a linear Fokker-Planck equation, the steady state probability distribution is known to be a Gaussian [3| . Thus we 
can describe the explicit form of P ss (r) simply as follows. 



P ss (r) 



Vdet[27rC(7)] 



exp 



(21) 



Here C(j) is the covariance matrix. 



C( 7 ) 



3 



{tqI? Tpj 
A(l + A) 1 + A 
T o7 



1 + A 




(22) 



For any physical quantity determined by R, the statistical average of a function of R in the steady state can be 
evaluated easily with the aid of the steady state probability distribution (f2"Tj) . Also, we can construct the constitutive 
equation of the anisotropic mobility model from the Fokker-Planck equation (|19p. We show the constitutive equation 
and compare it with some conventional constitutive equation models in Appendix [XJ 
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III. RESULTS 
A. Shear Viscosity 

In our model, the stress tensor is expressed as follows. 

a(R) = ^T RR (23) 

where is the number density of polymer chains. The steady state statistical average of eq (|23[) can be evaluated 
easily by using eqs (|2Tj) and ([22]) . The shear viscosity at the steady state is expressed as 

T/C-V) = (24) 
7 

where (• • • ) ss denotes the statistical average in the steady state (under shear). Finally we have the following expression 
for the steady state shear viscosity 77(7). 

R z l 1 + A 

Here rjo is the Newtonian viscosity in the limit of 7 — > 0, 770 = vakBTTo/2. Eq (|25[) is a monotonically decreasing 
function of j 2 . Thus we find that our model can reproduce the shear thinning behavior. 

Here we note that the shear thinning mechanism of our model is somehow similar to ones of the Giesekus model 
[2^ | or the Johnson-Segalman (JS) model 0, [47]|. Both the Giesekus model and the JS model are kind of linear 
dumbbell model which describes the dynamic behavior of viscoelastic fluids. (We briefly compare our model with 
the Giesekeus model or the JS model in Appendix [S]) In the Giesekus model, the anisotropic mobility comes from 
the hydrodynamic interaction from surrounding polymer chains. The anisotropic mobility tensor depends on the 
average conformation, which leads the effective acceleration of the relaxation and the shear thinning behavior. (The 
hydrodynamic interaction for a linear dumbbell model and the resulting shear thinning behavior have been extensively 
studied (4§|-|50|.) On the other hand, the characteristic feature of the JS model is the slippage effect, which allows 
a dumbbell to locally slip and effectively reduces shear deformation. In both models, a dumbbell apparently relaxes 
faster as we increase the shear rate, independent of their detail mechanisms. Our model behaves in a similar way. 
However, in our model, we did not explicitly consider the origin of the anisotropic mobility tensor (1111) . (In our model, 
the mobility tensor is designed based on the MD results, whereas in the conventional models it is usually designed from 
specific kinetic interactions.) We also note that nonlinear elasticity models (such as the finite extensibility nonlinear 
elasticity models (36|) can reproduce similar shear thinning behavior. It is difficult to identify the molecular level 
mechanism of shear thinning behavior only from the shear viscosity data. 

The first normal stress difference coefficient can be calculated in a similar way. 



(&xx — &yy)ss _ V^ B Tt q 

T 2 ~ 2 A(l + A) 



^{Af) = ^" xx . "yy /ss = "^£± '0 _ t (26) 



This also exhibits the thinning behavior. However, the dependences on the shear rate of 77(7) and ^1(7) are different. 
If we employ the power law type model for A (eq (|12[0 . at the high shear rate region we have 

77(7) tx j~ a (7 -> 00) (27) 
^1(7) cx 7~ 2q (7^00) (28) 

The shear viscosity and the first normal stress difference coefficient for the power law type model (with a = 9/11 (5lj |) 
are shown in FIG. [2] 



B. Rheo-Dielectric Response Function 



In dielectric measurements, we impose the time-dependent electric field in the y-direction (shear gradient direction), 
E y (t), and we measure the y-component of the electric flux density density. For a polymer chain which has a type- A 
dipole, the electric flux density (strictly speaking, the y-component of the electric flux density) D y can be expressed 
as follows [1|. 

D y (R) = EocEy + 4tw jlR y (29) 
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where is the effective dielectric constant due to the fast dynamics which is not resolved in our model, and /2 is the 
effective dipole intensity per unit backbone length of a polymer chain. At equilibrium, the Green-Kubo type formula 
(which is also referred as the Cole formula) gives 

v(t _ o _ 2,(1^(111^^ . (*»«>^) ea (30) 

Here, for simplicity we assumed that the correction factor for the internal electric field is unity. In equilibrium, our 
model (which reduces to eq ©) gives a single Debye type dielectric relaxation function. The result is 

<p(t) = Aeo-e^" (31) 

TO 

where we defined the dielectric intensity Aeo as Aeo = ^irvajj? R 2 /iksT . 

Under shear, generally we cannot use the Green-Kubo type response formula (|3"0)l . Although it is already shown 
that the Green-Kubo type formula can be used reasonably as a good approximation (23[, here we calculate the rheo- 
dielectric response function exactly to investigate the model properties precisely. Recently, Baiesi, Maes and Wynants 
[H IH derived a linear response formula in nonequilibrium states. Their formula does not involve any unclear 
approximations, and it can be applied to various nonequilibrium systems. Moreover, it is expressed as a simple form 
which enables us to evaluate linear response functions easily. We show a simple derivation of the Baiesi-Maes- Wynants 
formula in Appendix [Bj The Baiesi-Maes- Wynants formula gives the following form as the rheo-dielectric response 
function. 



,-,2 



(32) 



After straightforward calculations, finally the rheo-dielectric response function becomes a single Debye type decay 
function as follows. 

tp(j,t) = Aeo-e-V™ (33) 

TO 

Now we find that the rheo-dielectric response function (1331) is independent of shear rate 7, and thus it coincides with 
the equilibrium dielectric response function, eq pip . Then the rheo-dielectric intensity Ae(^) and the rheo-dielectric 
relaxation time r £ (7) simply become 

Ae 7 = Ae - 34 

6KB-L 

Ted) = To (35) 



Both are equivalent to equilibrium forms. This result is consistent with experimental data [201 - 4221 ] . Strictly speaking, 
the term "relaxation" should be used to describe "an approach" from a nonequilibrium state to an equilibrium state. 
But in this work, for convenience, we also call "an approach" to a perturbed state to a nonequilibrium steady state 
as "relaxation". 

Experimentally, it is convenient to use in the frequency domain expression (Fourier transform) of the rheo-dielectric 
expression, rather than the time domain expression (|35]) [J. The real and imaginary parts of Fourier transform of 
(|33|) become as £'(7,0;) = Ae /[1 + (wro) 2 ] and £"(7, u>) — Aequjtq/[1 + (wro) 2 ]. 

These results are rather trivial, because the y-component of the Langevin equation, eq (|15j) . is closed and contains 
only Ry(t) (no R x and R z contamination) and thus it is independent of the shear rate. If dynamics of x- and y- 
components are coupled, the rheo-dielectric function can be affected by the shear rate (as shown in Appendix |C|) . 
But we should be careful that even if there are no difference between the equilibrium dielectric response and the 
rheo-dielectric response under shear, the system is subjected to shear flow and thus not in equilibrium. Actually, as 
we showed in the previous subsection, the shear viscosity drastically decreased even when the rheo-dielectric response 
is not affected. 



C. Parallel and Perpendicular Moduli 



We consider the situation where a small, time dependent deformation is imposed to the constant-rate simple shear 
flow ©. Experimentally, so-called parallel or orthogonal superpositions are often utilized. In the case of the parallel 
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superposition, the shear rate 7 is modulated in time. This can be interpreted that we impose the perturbation velocity 
gradient tensor which has only the £?/-component to the system. The excess contribution for the xy-component of the 
stress tensor is then measured. On the other hand, in the case of the perpendicular superposition, the a;y-component 
is fixed to be constant and a small zy-component velocity gradient is imposed. The zy-component of the stress tensor 
is measured as the linear response. In the absent of the shear flow, these response functions coincide with the shear 
relaxation modulus G(t). The linear response theory gives 



G(t - 1') 



1 



cr xy (t) 



3k B T 



R x {t')R y (t') 



9v k B T 
i? 4 



(R x (t)Ry(t)R x (t )R y (t )} eq 



In our model, eq (j36]) reduces to a single Maxwell model. 

G{t) = G e~ 2t/T ° 



(36) 



(37) 



where we defined the characteristic modulus Go as Go = vaksT. 

Now we calculate the linear response of the stress tensor to the perturbation velocity gradient tensor, under steady 
shear. Unfortunately, the Baiesi-Maes-Wynants formula is limited for perturbations which can be expressed as 
perturbation potentials, and cannot be used in this case. Here we utilize the formula given by Seifert and Speck (3~3 | 
instead. The Seifert-Speck formula can be applied for non-potential type perturbations. We show a brief derivation of 
the formula in Appendix[Bj The Seifert-Speck formula gives the following expressions for the parallel and perpendicular 
moduli. 



G||(7,*-0 = ^(^( 



Gx(7,*"0 



knT 



R(t)) 



a zy (R(t)) 



Ry{t'] 
Ry(t' 



1 



1 



— +A sa (7) au un lRy{t ) 



dt' 



dR x {t') 



dRzW) dT{R{t')) 



After straightforward calculations, we find eqs and (|39p reduce to simple Maxwellian forms as follows 

G||(7,t) =G e- (1+ ^ /ro 
G±(-y,t) =G Q e- 2t ' T0 



(38) 
(39) 



(40) 
(41) 



We find that the parallel modulus f)40[) differs from the equilibrium shear relaxation modulus G(t) while the per- 
pendicular modulus (|4"Tj) is identical to G(t). Both Gii(7, t) and Gj_(7,i) approach Go at the short time limit. The 
parallel and perpendicular relaxation times under shear are given by 



T ||(7) = 



TO 



1 + A 



T di) = y 



(42) 
(43) 



The parallel relaxation time depends on the shear rate in the same way as the shear viscosity. We find that, in 
our model, the effect of the shear rate to the parallel relaxation modulus is observed only as the acceleration of the 
relaxation time. As we increase the shear rate, the parallel relaxation time decreases. This is in contrast to the 
rheo-dielectric function where we have no effect of the shear rate. Theoretically, the parallel modulus depends on the 
shear rate because it involves the correlation function of the x-component. 

As in the case of the rheo-dielectric response, it is convenient to use the Fourier transformed response functions in 
the frequency domain 4]. From eqs (|40[) and (|41[) . the real and imaginary parts (storage and loss) moduli become 
G\{i,u) = Go( W T||( 7 ))V[1 + (wt||(7)) 2 ], Gjf( 7 ,a/) = G ^T||(7)/[1 + ( W T||(7)) 2 ], G' ± (ff,u) = G ^r /2) 2 /[l + ^r /2) 2 ], 
and G" (7, uS) = Gq(ojtq/2)/\\. + {ujtq/2) 2 ]. We show parallel storage and loss moduli for several values of r c j with 
the power law model (eq (fT2")l ) in FIG. [3] 

Experimentally, both parallel and perpendicular relaxation times, rii and t±, decrease as the shear rate increases, 
and the decrease is more systematic for tii [l2j]. Our model predicts the decrease of the parallel relaxation time, 
which is consistent with the experimental data. On the other hand, the perpendicular relaxation time in our model 
is independent of the shear rate. This is because the dynamics of y- and z-components are not coupled to one of x- 
component. Thus we conclude that our model can reproduce the parallel modulus qualitatively but cannot reproduce 
the perpendicular modulus. (This is due to the oversimplification.) Although the parallel or perpendicular moduli 
have been analyzed or explained mainly on the basis of constitutive equation models (13 , 152l - l54j | , as far as we know, 
there is no analysis/explanation on the basis of the linear response theory. 
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IV. DISCUSSIONS 



A. Anisotropic Mobility Tensor 

Although our model is too simple to apply practical analyses, it still specifies some characteristic features of polymer 
dynamics under shear. We may comment that the anisotropic mobility tensor model proposed in this work can be 
understood as a simplified model studied in a previous theoretical work on the rheo-dielectric response [23j . In 
the previous work, we have studied the rheo-dielectric response function of a linearized Langevin equation for an 
entangled polymer. The anisotropic mobility tensor model is linear and has the same properties as the linearized 
Langevin equation model. (The expression of the rheo-dielectric response function is not affected significantly by 
the shear rate.) Our anisotropic mobility model reinforces the validity of the linearized Langevin equation model. 
Further, we expect that the linearized Langevin model will reproduce similar properties as the anisotropic mobility 
model (such as the acceleration of the parallel relaxation time). 

To model the dynamics of entangled polymers under fast shear, currently the convective constraint release (CCR) 
model [55J is widely employed. The CCR model claims that the effective relaxation time is accelerated under shear 
flow, due to the enhancement of the constraint release. Theories or simulations which take account the CCR effect 
achieved success to explain or reproduce rheological behaviour of entangled polymers under shear [55l - l60j . While the 
CCR model originally gives the expression for the relaxation time, here we interpret the modification of the relaxation 
time as the modification of the mobility (or the friction coefficient). Then the CCR model is interpreted as the shear 
rate dependent mobility model. The expression of the CCR mobility becomes as follows. 



A.ccr(7) 



1 . PCCR /r >r>\ 
Co k B T 



Co" + I^ 7(iW 



(44) 



where /3ccr is a positive constant of the order of unity and (...) represents the statistical (ensemble) average. In the 
steady state, (RR) can be replaced by the steady state average (RR) SS , which is determined self-consistently. 

It is obvious that the CCR mobility (l4"4")l is isotropic, and thus accelerates chain motion in all directions. As a 
result, linear response functions such as rheo-dielectric response function or the parallel modulus are affected by the 
shear flow. Although experimental data of parallel moduli can be reproduced well by the CCR model [l5l - [l7j . the 
rheo-dielectric response data fl8l I20l - l22j cannot be reproduced. In other words, the isotropic acceleration by the CCR 
is not consistent with the experimental data for the rheo-dielectric responses. As we have shown in the previous 
section, the anisotropic model can naturally overcome this difficulty. Therefore we consider that the CCR model 
needs to be modified to reproduce the anisotropic chain motion (or the anisotropic acceleration). For example, we 
may introduce a phenomenological anisotropic relaxation time tensor instead of a scalar relaxation time. Then the 
resulting model will reproduce the rheo-dielectric function which is insensitive to the shear rate, as well as the shear 
thinning. One simple possible modification is shown in Appendix [D] 

To compare the CCR model with our anisotropic mobility model from a different aspect, here we consider a general 
form for the mobility tensor. We consider the mobility tensor under shear flow in general flow and gradient directions. 
From the symmetry of the Langevin equation under rotational transform, the mobility tensor A should be invariant 
for rotational transform. That means, we can expand A into a power series of scalar and tensor invariants if the shear 
rate is not high. Then, up to the second order in 7, the expansion form of A is given as 

A(k) = -!- 1 + Li(k + «*) + L 2 tr(K • + L 3 k • k* + L 4 k' ■ k + 0(j 3 ) (45) 
Co L J 

Here {Li} is a set of expansion coefficients. We find that the CCR mobility model (|44|) can be reproduced by setting 
L\ = L3 = L4 = 0, while our anisotropic mobility model can be reproduced by setting L\ = Li = L4 = 0. (If 7 is 
sufficiently small, we can set (R x R y ) = 7T0-R 2 /6 + 0(j 2 ) in eq (|44j) . If k is given by eq ©, only the xcc-component 
of k • is nonzero.) Therefore both the CCR model and our model are allowed from the symmetry argument. 
The mobility tensor should be modelled so that the resulting dynamics reproduces required properties (such as the 
insensitivity of the rheo-dielectric response function to the shear rate). We may employ a different mobility tensor 
model which is reduce to eq (|45|) . For example, by setting L2 — Lj, = L4 = 0, we have a nondiagonal mobility 
tensor model which takes account of the kinetic coupling effect. (See Appendix 0) Here we should note that in 
the conventional approach [30|, the mobility tensor is expressed as a function of the average conformation tensor 
(not as a function of the velocity gradient tensor), and thus it is expanded into a power series of the conformation 
tensor. However, from the view point of the nonequilibrium statistical physics [3lj . in principle, the mobility tensor 
can depend on the velocity gradient tensor and can be modelled as eq (fTTj) or eq f|45|) . (We should also note that the 
current approach is limited for simple shear flows, and for other flows, such as elongational flows, we need to construct 
the mobility tensor model.) 
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One may argue that if an isotropic mobility (which can be utilized near equilibrium) should be replaced by an 
anisotropic mobility for nonequilibrium systems, then an equilibrium free energy (J3]) should also be replaced by a 
nonequilibrium effective free energy. This is indeed correct in general. Nevertheless, for polymeric systems we can use 
the equilibrium free energy §3§ safely even under fast shear. This is because in polymeric systems, the stress-optical 
rule [j, |6l| is known to be hold even under fast shear. (Although it is known that the stress-optical rule fails for some 
cases, such as systems under fast extensional flow 6^, 63[, it is valid under the current situation.) The stress-optical 
rule relates the chain conformation and the force exerted by the chain, and holds only when the force is proportional 
to the bond vector. This gives the linear entropic elasticity model for the free energy, and as long as the stress-optical 
rule holds we can justify the use of the equilibrium free energy ((3|) even under fast shear. 

At the end of this subsection, we may comment on one assumption used in our model. We have assumed the 
fluctuation-dissipation like relation for the thermal noise (eq (jlOp ). Such a relation does not necessarily hold in the 
nonequilibrium states, and thus we can employ nondiagonal mobility tensor model. From the view point of the linear 
response theory, whether the mobility tensor is diagonal or not is not so essential (see Appendix [Bj. Although we do 
not discuss further in detail about the nondiagonal mobility tensor because it is beyond the scope of this work, we 
expect it may be required to describe polymer dynamics precisely under general flow conditions. 



B. Entropy Production and Steady State Probability Current 

We have derived expressions for several linear response functions. In several pieces of theoretical work, the violation 
of the fluctuation-dissipation relation is interpreted as the entropy production rate [64], [65| , the steady state probability 
current 66 68], or other related physical quantities. Here we calculate the entropy production rate for our model and 
investigate how it is related to the response functions. 

Following the standard definition [69i - l71| we define the steady state entropy production rate per unit volume as 
follows. 



S(7) 



T 



dr 



J.(r).A- 1 (7).J»(r) 



P ss (r) 



where J ss (r) is the steady state probability current defined as 

~dT{r) 



J ss (r) = -A( 7 ) 



dr 



-P ss (r) + k B T 



dP ss (r) 



Or 



k ■ rP ss (r) 



(46) 



(47) 



The steady state probability current (|4"T)) satisfies the steady state condition, (d/dr) ■ J ss (r) = 0. After the straight- 
forward calculation, we have the explicit expression for the steady state entropy production rate in our model. 



S(7) = ^A(7): J drP ss (r) 



3k B T 



K 2 



r — A ■ k ■ r 



3k B T 



r — A 1 (7) • k ■ r 



3^ fc| r 

R 2 



Mi) ■■ i 



"oCgA ( 1 
T 



T Cxx 



27 7 2 
Ln _l —C 

~ <~^xy T ~ ^yy 
T A X 2 



Tt 2 



(48) 



to A(l + A) 



As expected, the entropy production rate (|48j) is the even function of 7, and it is nonzero unless 7 = 0. We show the 
entropy production rate for the power law type model (a = 9/11 and t c = to) in FIG. HI We can employ another 
definition for the entropy production rate, S = A f{<J X y)ss/T [72]. This gives a slightly different form from eq (|48p. but 
the result is qualitatively the same. 

Some nonequilibrium linear response theories state that, the entropy production rate (which may be interpreted as 
the distance from equilibrium) is related to the violation of the Green-Kubo type response formulae. We call such 
a picture as the entropy production picture. In the entropy production picture, one expects that if the system is 
not in equilibrium, there should be the correction terms in the linear response formulae which is directly related to 
the entropy production rate. However, as we have already shown, the rheo-dielectric response function is unchanged 
even under shear in our model. Besides, in the previous work [23], the Green-Kubo type formula was shown to be 
approximately valid for the rheo-dielectric response. These results mean that, some linear response functions do not 
change their forms even in the nonequilibrium states. This may sound inconsistent with the entropy production 
picture. This is because what appears in linear response formulae is not the entropy production rate itself but the 
derivative of the entropy production rate with respect to an external perturbation field. 
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The situation may be clearer if we employ another picture, which focus on the steady state probability current. 
A linear response function in nonequilibrium steady state can be expressed as the sum of the Green-Kubo type 
equilibrium form and the correction term, which involves the probability current. We call this picture as the Lagrangian 
moving frame picture (6(1 . The Lagrangian moving frame picture gives the following response formula. 



K AB {t-t') 



1 



kgT dt 



1 



k R T 



A(t)v m (R(t')) 



dR(t') 



(49) 



where v ss (r) is the mean steady state streaming velocity defined as v ss (r) = J ss (r) / P ss (r) . The second term in the 
right hand side of cq (|49"1) is the correction term due to the nonzero steady state probability current. The correction 
term can be zero even if «7 ss (t") ^ 0. In the case of the rheo-dielectric response in our model, the correction term is 
exactly equal to zero. This can be shown straightforwardly. 



Ry(t) 



R x (t>) - -P^Ry{t') 



= e ~(t-t')/T 



= 



1 + A 

{RyRx}sS ~ 



TQ 7 

1 + A 



(Ry) S 



(50) 



Moreover, it is quite difficult to separate an experimentally measured response function into two terms as eq (|49l) . 
This is because the first term in the right hand side of eq (|49|) (the Green-Kubo type term) can also depend on the 
shear rate. 

From the results and discussions above, we consider that even if we measure the rheo-dielectric function and the 
entropy production rate of the same system under shear simultaneously, we will not be able to verify the violation of 
the Green-Kubo type formula precisely. Thus we consider that the entropy production or the Lagrangian moving frame 
pictures are not so useful to analyze rheo-dielectric responses or other linear responses under shear. To investigate 
dynamics of polymer chains under shear in detail, we consider it is better to measure other linear response functions, 
such as the parallel modulus, instead of the entropy production rate (or the corresponding heat flow). Combination 
of several linear response functions will provide us detail information about the dynamics of polymer chain [5-9]. 

To be fair, we should mention that for microscopic systems (such as a colloid particle driven by an optical trap 
[zHlzIl]) the entropy production rate or related quantities can be utilized successfully to characterize the nonequilib- 
rium features. For carefully designed microscopic systems we can measure correlation functions or the steady state 
probability current directly by microscopes. These physical quantities are essential in the entropy production or La- 
grangian moving frame pictures. But in macroscopic systems, it is quite difficult or practically impossible to measure 
several physical quantities. Thus a different approach for macroscopic systems is naturally required. 



C. Dependence on Architecture of Polymers 



We have shown that our anisotropic mobility tensor model can reproduce rheo-dielectric response behavior quali- 
tatively. That is, the rheo-dielectric response functions of linear polymers are insensitive to the shear rate. However, 
rheo-dielectric response functions of star polymers are reported to slightly depend on the shear rate [2(| ■ This cannot 
be explained by our model. In this subsection, we consider why our model fails to describe star polymers and possible 
ways to improve the model. 

Arsac et al [28[ fit experimental rheology data to the JS model and determined the slip factors (the fitting parameters 
in the JS model). They found that the slip factors for linear polymers are nearly independent of the flow regime 
(transient or steady state) or the molecular weight distribution. This means that the JS model can reproduce 
dynamics of linear polymers in spite of its very simple form. However, for branched polymers the slip parameters 
depend on various factors. We may say that the dynamics of entangled linear polymers is rather simple, in a sense. 
Thus we consider that the dynamics of star polymers cannot be described well by a simple model like the JS model. 

Matsumiya, Watanabe and coworkers [6-9] measured and analyzed the linear viscoelasticities and dielectric re- 
sponses of entangled linear and star polymers. They quantitatively tested the dynamic tube dilation (DTD) model 
[75j for linear and star polymers. They reported that the simple DTD model explains the experimental data for 
linear polymers well @, while it fails for star polymers 043- This failure is attributed to the overestimate of the 
equilibration by the constraint release at the long time region [76T - |78| , or the events that newly created entanglements 
push out the old entanglements toward chain ends [791 ]. These experimental results indicate that a rather simple 
model can describe the dynamic behavior of entangled linear polymers but not of star polymers. 
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Thus we expect linear polymers have rather simple dynamic properties even under fast flow while star polymers 
do not. This implies that our simple anisotropic mobility model will be able to describe dynamics of linear poly- 
mers qualitatively well. On the other hand, the dynamics of branched polymers (including star polymers) is much 
more complicated compared with linear polymers and too simplified models (like ours) cannot explain dynamics of 
branched polymers. Then we can conclude that our model should not be applied for star polymers directly, and some 
modifications are required. 

There are several possible ways to improve our model. For example, we can employ the nondiagonal mobility tensor 
model, which represents the kinetic coupling between dynamics of different directions (the kinetic coupling model). 
Such a model can reproduce the dependence of the rheo-dielectric response to the shear rate to some extent. We 
show the rheo-dielectric response function for a simple and weak kinetic coupling model in Appendix [Cj Similarly, we 
can employ the conformation dependent mobility tensor model |25l - [27l |8CH82| • The conformation dependent mobility 
kinetically couples the dynamics for different directions, and will give similar results as simple kinetic coupling models. 

Another possible way is to employ a fine scale description such as the full bead-spring type model with topological 
constraints (83[. Integrating our anisotropic mobility model into bead-spring type models will allow us to study the 
dependence of the rheo-dielectric behavior on the polymer architecture. Anyway, the rheo-dielectric behavior and 
dynamics of entangled star polymers are still not fully understood and further theoretical developments are required. 

V. CONCLUSION 

In this work, we proposed the anisotropic mobility model for polymers under shear. In the anisotropic mobility 
model, the mobility tensor or the friction coefficient tensor becomes anisotropic and dependent on the shear rate. The 
anisotropic mobility tensor model is consistent with NEMD simulation results and thus we expect that our model 
captures qualitative and essential nature of polymer dynamics under shear. 

We calculated the shear viscosity, the rheo-dielectric response function, or the parallel and perpendicular moduli. 
Our model gives the rheo-dielectric function which is independent of the shear rate, even when the shear rate is 
sufficiently high and shear thinning is exhibited. This is qualitatively consistent with the experimental results. Our 
model gives the parallel relaxation time which decreases with increasing the shear rate. This is also qualitatively con- 
sistent with the experimental results. Of course, the shear-rate insensitive rheo-dielectric relaxation and acceleration 
of the parallel relaxation observed in experiments may result from not only the anisotropic mobility but also from 
other factors (such as full DTD in the linear response regime). However, the current study demonstrates that the 
anisotropic mobility could play an important role in the relaxation processes. 

To examine the properties of our model in detail, we compared our model with other models or theories. We compare 
our model with the CCR model. Both our model and the CCR model accelerate the dynamics of polymers under 
shear. Our model accelerate the dynamics anisotropically while the CCR model accelerate the dynamics isotropically. 
Judging from the experimental results, we consider our model is more suitable to describe the dynamics of polymers 
under shear. 

We also compared our result with the recent linear response theories for nonequilibrium systems. Although the 
entropy production rate or the steady state probability current are widely utilized in recent models, we showed that 
they are not so useful to analyze or understand the rheo-dielectric response function. We consider that the combination 
of several different linear response functions will be reasonable to investigate polymer dynamics under shear. 

Although our model can explain the essential feature of linear polymers under shear, it should be improved or 
modified further. For example, our model can not explain the experimental results for star polymers under shear. We 
did not explicitly consider the effect of entanglements, which will be important for star polymers. The integration of 
our anisotropic mobility model into fine scale models or the generalization of our model to general flow conditions is 
considered to be an interesting subject of future work. 
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Appendix A: Constitutive Equation for Anisotropic Mobility Model 

In this appendix, we derive the constitutive equation from the anisotropic mobility model and compare it with some 
conventional models. For simplicity, we assume that the system is homogeneous in this appendix. (The generalization 
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for inhomogeneous systems is straightforward.) We consider to express the constitutive equation as the dynamic 
equation for the following time-dependent conformation tensor. 



C(t) = (R{t)R(t)) = J drrrP(r,t) 



(Al) 



The time evolution equation for the conformation tensor can be easily calculated from the Fokker-Planck equation 
(|19p . After a straightforward calculation, we have the following constitutive equation for the anisotropic mobility 
model (eq |(HJ) together with eq (|TT|) or eq (1451) ). 



v 

C(t) 



3k B T 
~R2- 



[C ■ A(k) + A(k) • C] + 2k B TA(K) 



(A2) 



where we defined the upper-convected derivative as C = dC / 'dt— n-C ' C '• k 1 . It should be noticed that the anisotropic 
mobility model is designed around the steady state under simple shear and thus the corresponding constitutive equation 
(|A2[) is also applicable around the steady state. (It is not suitable to calculate, for example, the start-up shear flow. To 
study such transient phenomena, we will need to describe the time evolution of the the mobility tensor A explicitly.) 

It is informative to compare eq ()A2j) with other constitutive equation models. Although there are many constitutive 
equation models for polymeric systems [30j . for the sake of simplicity, here we limit ourselves to rather simple models. 
One of the simplest constitutive equation models with anisotropic mobilities is the Giesekus model. The Giesekus 
model (27j employs the conformation tensor dependent mobility, whereas the anisotropic mobility model employs the 
mobility which does not depend on the conformation tensor. The Giesekus constitutive equation can be expressed as 
follows. 



v 

C(t) 



Co 



(1-0)1 + 



3a 



3k B T 



C + k B Tl 



(A3) 



where a is a phenomenological constant (0 < a < 1). (The Giesekus model corresponds to the pre- averaged and 
linearized Curtiss-Bird model [25l - [27| .) Eq (IA3|) can be obtained by replacing A(k) in eq (|A2I) by [(1 — a)l + 
(3a/R 2 )C]/( Q . Inversely we can replace [(1 - a)l + (3a/R 2 )C}/( by A(k) to obtain eq (|X2]l from eq (|A3]l . One 
may interpret the anisotropic mobility model as the conventional anisotropic tensor model with some approximations, 
for example, the pre-averaging around the steady-state. (We note that it is not simple to analyze eq (|A3|) due to its 
nonlinearity, unlike eq (|A2I) . The conformation tensor independent mobility tensor makes analyses in the main text 
simple and tractable.) 

Another simple constitutive equation model is the Johnson-Segalman (JS) model |46[. The JS model employs the 
Gordon-Schowalter derivative |84j to produce the non-affine motion, which can be interpreted as the slippage effect. 
Here we rewrite the JS model by using the upper-convected derivative, to compare it with eq (|A2I) . 



v 2 
C(t) = ± 
4o 



3k B T 



C + k B Tl 



1 



[(k + k*) -C + C-(k + K*)] 



(A4) 



Here a is a phenomenological constant (— 1 < a < 1), which is sometimes called the slip parameter. We find that the 
form of the JS model ()A4j) is somehow similar to the anisotropic mobility model (|A2|) . The last term in the right 
hand side of eq (|A4[) is mathematically similar to the first term in the right hand side of eq (|A2[) . Thus we expect 
that the anisotropic mobility model and the JS model will show qualitatively similar dynamical behavior in some 
cases. However, the origin of that term in the JS model is the non-affine motion (or the slippage). Unlike the case of 
the Giesekus model, we cannot obtain the anisotropic mobility model (IA2[) by simply replacing a part (such as the 
mobility tensor) in eq (|A4I) . 



Appendix B: Derivation of the Baiesi-Maes-Wynants Formula by the Path Integral Formalism 

In this appendix, we show the derivation of the Baiesi-Maes-Wynants formula [H, H3| in steady state based on the 
path integral formalism. The derivation of the linear response formula in nonequilibrium steady state based on the 
path integral formalism is first shown by Seifert and Speck [34l |. Here we mainly follow their derivation. It is worth 
noting that the path integral formalism had already utilized by Ohta and Ohkuma [68[ to derive a similar but slightly 
different formula, in prior to the Seifert-Speck theory. 

The Seifert-Speck formula reduces to the Baiesi-Maes-Wynants formula if the perturbation is expressed as an 
perturbation potential. As far as we know, an explicit derivation of the Baiesi-Maes-Wynants formula by the path 
integral formalism has not been presented. 
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Here we consider a general nonequilibrium system, of which dynamics follows the Langevin equation. We denote 
the dynamic variables which obey the Langevin equation as Xi , A" 2 , ■ ■ . , X n (with n being the number of independent 
stochastic variables). For convenience, we use the Ito stochastic calculus [85| for the stochastic differential equation. 
(One can employ the Stratonovich calculus instead of the Ito calculus. Although the calculation below becomes 
somehow complicated, the result is essentially the same.) For simplicity we introduce the shorthand notation X = 
(Xi, X2, • • • , X n ). Because many physical quantities depend on X(t), we also introduce a shorthand notation for a 
function of X(t) as f(t) = f(X(t)). The Langevin equation for X(t) is described as 

dX(t) w , „ a 



dt 



V(t,h(t)) + k B T 



dX(t) 



(Bl) 



where V(t, h(t)) is the average change rate of X (which may be understood as a sort of velocity), h(t) is the external 
perturbation, A is a positive definite symmetric tensor, and £(t) is the Gaussian noise. V(t, h{t)) is decomposed into 
the reference part which is independent of h(t) and the perturbation part which is linear in h{t). 



V(t,h(t)) = V (t) + V 1 (t)h(t) 



(B2) 



where Vo(t) is the average velocity exerted by the interaction potential or the external force, and V\(t)h(t) is the 
perturbation term. We assume that h(t) is sufficiently small to neglect higher order terms in h(t). The third term 
in the right hand side of cq (IB1[) is the stochastic drift term which cancels unphysical probability current £(i) 
satisfies the following equations. 



(€(*)> = 

<£(*)£(*')> = 2k B TA(t)S(t - 



(B3) 
(B4) 



Eqs (|B3I) and (|B4I) can be interpreted as the fluctuation dissipation type relation. Or, inversely we can define k(t) 
via eq (|B4l) . 

The probability that a trajectory X(t) is realized, which we may call the path probability (or the path weight), 
can be calculated from the distribution of the noise £(t). Because £(t) obeys the Gaussian distribution, the path 
probability 7 :> [X(-)] can be calculated as follows (86| . 

dX{t) V{t,h{t))-k B T-^—.k{t) 



V[X{-)\DX =Afexp 
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dX(t) 
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dX{t) 



k{t) 



2k B T 
dX(t) 



dtVx{t) ■ A-\t) 



dt 



V {t)-k B T- 



d 



k(t) 



h(t) + 0{h 2 
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V [X(-)]VX 



(B5) 



8X{t) 

where Af is the normalization factor and "Po[-^(')] ^ s ^ ne P a tb probability at the reference state (without any pertur- 
bations). The time derivative is interpreted as the retarded derivative [86| (which is consistent with the Ito calculus), 
and thus JV is independent of h(t). The path probability at the reference is defined as 

~dX{t) w,., , „ d 



V [X(-)]VX=J\fexp 
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4fc s T 
dX{t) 
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dX(t) 



A(t) 
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(B6) 



VX 



dt ' uv " y ox(t) 

By using the path probability (|B6|) . we can define the steady state statistical average as the following path integral. 



<->. 



J VX---T [X(-)} 



(B7) 



Since we are interested in the linear response, the 0{h 2 ) term in eq (|B5|) can be safely neglected. Then, the 
statistical average of a physical quantity A at time t with perturbation can be expressed as 



(i(t))= VX A{t)V[X{.)} 



2k B T 



dt' ( A(t)Vi(t') ■ A _1 (t') 



dX(t<) * d 
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(B8) 
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Here, A ss = (A(t)) ss is the steady state statistical average at the reference state (without perturbation). From time 
translational symmetry, A sa is independent of time t. We have the following expression as the response function of 
A(t) to the perturbation h(t') from eq (|B8|) . 



Tl A (t - t 1 ) 



1 



2k B T 



dX(t' 
dt' 



V (t')-k B T 
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A(f) 



(B9) 



Especially, if the perturbation is caused by a perturbation potential, the perturbation term can be rewritten by 
using just a single scalar quantity. If we assume the fluctuation-dissipation type relation, then Vi(t) can be rewritten 
as follows. 



V 1 (t)=A(t) 



dB{t) 



dX(t) 

Where B is the scalar quantity which is conjugate to h. Eq (|B9[) can be simplified as follows. 
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where we have utilized the Ito formula 

dB(t) 

dt ~ dt ' dX(t) ' " v " ; ' dX{t)dX{t) 
and defined the backward generator £' (which has the same form as the associate Fokker-Planck operator) as follows. 
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Eq (|B11| is nothing but the Baiesi-Maes-Wynants formula [32', 33]. Although eq (|B11|) is simpler than eq (|B9j) . we 
should notice that eq (jBllj) can be utilized only when the perturbation is given as a perturbation potential and the 
fluctuation-dissipation type relation holds. We should directly use eq (|B9[) if these conditions are not satisfied. 



Appendix C: Weak Kinetic Coupling Between Different Directions 



In this appendix, we consider the kinetic coupling effect between the x- and y-direction dynamics. From the NEMD 
simulation results |4Cl |41| , we expect that the iry-element of the mobility tensor is sufficiently small compared with 
the diagonal elements. We can employ the following nondiagonal mobility tensor model as a simple kinetic coupling 
model. 



A(7) = 



1 aro7 
ar 7 1 







1 



(CI) 



Here a <c 1 is a parameter which represents the coupling strength. (Eq (|C1|) is obtained by setting L2 = L3 = L4 = 
in eq (I45[l .) Unlike the model considered in the main text, this mobility model does not accelerate the dynamics 
in x-direction explicitly. Nonetheless, this model can be used to demonstrate how the kinetic coupling affects the 
rheo-dielectric behavior. 

We consider a as the perturbation parameter, and expand physical quantities into the power series of a. To see 
how the kinetic coupling affects the viscoelastic or dielectric properties, it is sufficient to consider only the leading 
order terms (which are proportional to a). First we consider the steady state probability distribution. Since the 
Fokker-Planck equation is linear, the steady state probability distribution can be expressed as a Gaussian. Then the 
covariance matrix can be expressed as follows. 
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(R y (t)R x (0)) ss = C^e-*/™ + O(a) 
From the Baiesi-Maes-Wynants formula, the rheo-dielectric response function becomes as follows. 
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The rheo-dielectric intensity can be calculated to be 

POO 

A £ ( 7 )= / dt<p{t,j) = Ae [l- a (ro7) 2 + 0(a 2 )] 



(C4) 
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(C6) 



From eq (|C6p we find that the rheo-dielectric intensity is decreased by the kinetic coupling. This is in contrast to 
the case of the anisotropic mobility model, where the rheo-dielectric intensity is independent of the shear rate. By 
performing the Fourier transform for cq (|C5[) . we have the following real and imaginary parts of the rheo-dielectric 
response function in the frequency domain. 
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From eqs (|C7|I and (|C8|) . we find that the effects of the kinetic coupling and the shear rate can be represented by an 
effective coupling constant, a(r 7) 2 . We show e'(j,v) and e"(j,uj) by eqs (|C7I) and (|C8|) in FIG. [5j We can observe 
that the dielectric loss e"( 7 , oS) for tow < 1 decreases as the effective coupling constant a(ro 7 ) 2 increases. On the 
other hand, we can observe that the rheo-dielectric response functions are insensitive to the shear rate for tol o > 1. 
These are qualitatively similar to experimentally observed rheo-dielectric behavior of entangled star polymers [2(| . 



Appendix D: Anisotropic Version of the Convective Constraint Release Model 

As we discussed in the main text, the convective constraint release (CCR) model isotropically accelerates the 
dynamics of a polymer chain. In order to reproduce the rheo-dielectric response behavior, we need to modify the 
CCR model to be anisotropic. In this appendix, we consider a simple anisotropic version of the CCR model, as a 
possible modification for the CCR model. 

The key feature of the CCR model is that the characteristic relaxation time is accelerated by the product of the 
velocity gradient tensor and the chain conformation tensor. To make the CCR model anisotropic, here we employ the 
following mobility tensor. 

Accr(7) = y X + IN [« ' « RR ^ ( RR ^ + « RR ^ ( RR ^ ■ K ^ ( D1 ) 

where /?ccr ' s a positive constant of the order of unity and the velocity gradient tensor k is given by eq ([6]). Eq 
(jDljl is similar to eq (|44|) but the scalar diadic product k : (RR) is replaced by second rank tensor products. For 
the steady state under shear, the ensemble average (RR) is replaced by the steady state average {RR) SS , and we can 
obtain the steady state probability distribution explicitly. The anisotropic CCR model (|D1[) has the Gaussian form 
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steady state probability distribution ([21)) with the following covariance matrix. 



It is straightforward to show that eqs 
viscosity is expressed as follows. 
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21]) . (|D1|) . and (|D2j) satisfy the steady state condition (|20|) . The steady state 



2[V^cr(^o7) 2 + 1-1] 
/W( r o7) 2 



(D3) 



Eq (|D3I) is a monotonically decreasing function of 7 2 , and thus the anisotropic CCR model shows the shear thinning 
behavior. Eq (IL)3[) is similar to that for the steady state shear viscosity in the original CCR model. At the high shear 
rate region, eq (|D3|) approaches to the following asymptotic form. 



nil) -+ vo 



(7 — > 00) 



As in the original CCR model (55[, eq (|D4[) becomes consistent with the Cox-Merz rule when we set P[ 



CCR 



(D4) 



4. We 



can straightforwardly show that the first normal stress difference coefficient by the anisotropic CCR model is also 
similar to one by the original CCR model. 

The steady state covariance matrix (|D2I) looks similar to the steady state covariance matrix for the anisotropic 
mobility model (122[) . In fact, the anisotropic CCR model reduces to the anisotropic mobility tensor model 
Substituting eq (ID2[) into eq (|D1I) . we have 
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(D5) 



This has the same form as the power law type anisotropic mobility tensor model (eqs (jlip and ([12"]) ') with r c = ■J (3' ccr tq 
and a = 1. Then it is clear that the anisotropic CCR model successfully reproduces essential properties of our 
anisotropic mobility tensor model. However, it should be noticed that the anisotropic CCR mobility tensor (|D1[) 
depends on (RR) in the presence of the external perturbation field. Because (RR) generally depends on the applied 
perturbation field implicitly, the anisotropic CCR mobility tensor gives additional contributions for linear response 
functions. Thus the linear response properties of the anisotropic CCR model will be slightly different from ones of 
our anisotropic mobility model. We should carefully calculate the contribution from the mobility tensor to get the 
linear response functions for the anisotropic CCR model. (Fortunately, even if the mobility tensor depends on the 
perturbation field, the path integral formulation shown in Appendix IBI is still valid.) 
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FIGURE CAPTIONS 

FIG. [T] Schematic image of a polymer chain under shear. The solid curve and the solid arrow represents the polymer 
chain and its end-to-end vector of the chain, respectively. The dotted arrows represent the flow field given by eq ([5]) ■ 
R and 7 are the end-to-end vector and the shear rate, respectively. 

FIG. [2] The shear viscosity 77(7) and the first normal stress difference coefficient calculated by the power law 

type model (eq (|12"j) ). The exponent is set to a = 9/11. 

FIG. [3] The parallel storage and loss moduli GUjfUj) and G',!(j,cj), calculated by the power law type mobility model 
(eq (fl2)) ). The exponent is set to a = 9/11. 

FIG. S]The entropy production rate £(7) calculated by the power law type model (eq (|12p). The exponent and the 
characteristic crossover time are set to a = 9/11 and t c = to, respectively. 

FIG. [5] The rheo-dielectric response functions for the weak kinetic coupling model. The effective coupling constant is 
varied as a (to 7) 2 = 0, 0.1, 0.2. The case of a (to 7) 2 =0 corresponds to the equilibrium dielectric response. 
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